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Abstract 

The results from developing and applying the notions of functional self-similarity 
and the Bogoliubov renormalization group to boundary-value problems in mathe- 
matical physics during the last decade are reviewed. The main achievement is the 
regular algorithm for finding renormalization group-type symmetries using the con- 
temporary theory of Lie groups of transformations. 

1 Introduction 

The notion of functional self-similarity (FS) was introduced in mathematical physics by 
one of the authors in the early 1980s [Q] (see also P, 13]). The basis of this introduction 
is that solutions of a wide class of problems analyzed by the renormalization group (RG) 
method are invariant w.r.t. the group transformations that involve not only natural in- 
dependent variables of a problem but also parameters of boundary conditions imposed 
at some "reference" point. The RG transformation then corresponds to reparameteriz- 
ing a solution by changing (shifting or rescaling) an independent (coordinate) variable 
while simultaneously performing a functional transformation of (boundary) characteris- 
tics of the selected functions when passing to another reference point. The corresponding 
"transformation functions" are governed by group functional equations. 

More precisely, we consider the so-called renormalization transformations (the Dyson 
transformations) in quantum field theory (QFT), which constitute a continuous one- 
parameter group, i.e., the Lie group of transformations (if we use a standard mathematical 
language). This group was discovered |4| and used [|J to analyze QFT singularities. We 
call this group the QFT RG or the Bogoliubov RG to distinguish it from the approxi- 
mate RG, which was introduced by Wilson |]] to analyze critical phenomena in statistical 
physics problems. 
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The reasoning above pertains to the RG transformations considered within the QFT 
approach. In the middle 1980s, the RG method became widely applied to classical prob- 
lems in mathematical physics; this was initially connected with transforming and applying 
methods developed for QFT and statistical physics. Such applications were based, first, 
on the fact that physical problems manifest the FS property, which permits segregating 
the characteristic variables of a problem (independent and dynamic variables, parameters, 
and boundary data) and finding transformations that preserve the solution, and, second, 
on an advantageous mathematical description that uses only differential and integral- 
differential equations to describe the model. The RG approach is also favorable because 
the corresponding RG-type transformations can be constructed using regular methods, 
which are used to find symmetries in the group analysis of differential equations (DE). 
Substantial progress has been achieved in group analysis since the RG method first ap- 
peared in theoretical physics several decades ago. 

The transition to boundary- value problems in mathematical physics enriched the 
mathematical content while preserving the main property of the FS, i.e., the solution 
of a physical problem is invariant w.r.t. a special class of transformations. The initial 
notion of the FS transformations as point transformations of independent variables was 
recently generalized to contact transformations, transformations defined by formal series, 
etc. Because the infinitesimal transformation approach is convenient for the group anal- 
ysis of mathematical models based on DEs, we can also formulate the FS in terms of 
infinitesimal operators determining the corresponding RG symmetry (RGS). The oper- 
ators then appear as a result of the standard RGS construction procedure, and the FS 
condition arises in the course of this procedure. 

We review the evolution of the FS notion connected with the implementation of the 
RGS method in mathematical physics and describe the results obtained in the RGS frame- 
work. In Sec. 2, we introduce both infinitesimal and finite RG transformations, establish 
their connection with the mathematical physics notion of powerlike self-similarity (auto- 
modelness), and introduce the notion of the FS. We then construct RGSs in mathematical 
physics. The FS property appears as an ingredient of this construction, and we trace how 
the FS notion changes when passing from the Bogoliubov RG to mathematical physics 
models. 

In Sec. 3, we consider examples of the FS transformations that are close to the QFT 
transformations because they are realized as groups of one-parameter transformations 
w.r.t. some independent variables. We present examples where the group extends and 
generates FS transformations with operators constituting a finite-dimensional algebra. 
Because an ODE or a system of such equations governing the mathematical model dis- 
cussed in this section coincides with the Lie equations arising in the QFT applications of 
the RG method, the formalism developed for ODEs can also be applied to Lie equations. 

In Sec. 4, we collect examples of FS conditions written either as a first-order PDE or as 
higher-order differential relations. We discuss a novel form of the FS transformations that 
can be presented as formal infinite series, not as traditional algebraic relations. The FS 
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transformations acquire such a form when the RGS are of the Lie-Backlund group type,[] 
and the FS condition claims that the particular solution of a boundary-value problem is 
invariant w.r.t. these transformations. 

In Sec. 5, we analyze the FS conditions for systems with small parameters. The 
FS conditions in such systems permit finding approximate RGSs and using them to solve 
problems with arbitrary boundary conditions. We present examples of approximate RGSs 
for several problems in nonlinear physics. 

In Sec. 6, we describe the perspectives of the FS and its use to construct RGSs in a 
broad class of mathematical physics problems. 



2 The QFT RG and the FS condition 

2.1 Simple Bogoliubov RG transformations 

To illustrate the FS notion, we consider the simplest RG transformation, which is a 
simultaneous one-parameter transformation of one independent variable in a problem 
(the coordinate £, for example) and a characteristic g of a solution, 

T{\):{i^£' = i-\, g^g> = G{\,g)} , G(0,g) = g. (1) 

The function G(£,g) satisfies the functional relation 

G(£ + X,g) = G(£,G(X,g)), (2) 

which corresponds to the group composition law for the transformation operators, T(Ai) • 
T(A 2 ) =T(Ai + A 2 ). 

Transformation operator fll]) can sometimes be conveniently represented explicitly by 
writing a transformation of a function F(£, g) in the form 

T(\)F(£,g) = e- XR F(£,g) (3) 

using the infinitesimal RGS operator (or the RG operator) 

R = d i -P(g)d g , p(g) = ^^l , (4) 

OA A=0 

whose coordinate f3(g) is the derivative of the function G at A = 0. The infinitesimal RG 
transformations can be conveniently represented via the operator R; the group nature of 
the operator T(A) = e~ XR is obvious from "finite shift" representation (|3]). 
The condition 

RI(£,g) = d e I - (3{g)d g I = , (5) 

lr The set of group variables of the generalized Lie group of transformations (the Lie-Backlund group) 
includes derivatives of the desired functions w.r.t. independent variables. 
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determines an invariant of the RG transformation, which is the function I(£, g) = I(G(l, g)) 
= I (0,G(£, g)) of one argument in this case because the transformation function G(£,g) 
is an invariant itself by virtue of Eq. (|2[). The condition 

RC(£,g) = <p{g)C(e,g) (6) 

determines a covariant, i.e., a quantity that transforms according to a representation of 
the FS group or the RG. Covariants are important for QFT applications of the RG. 

The representation of a RG transformation via an infinitesimal operator R is equiva- 
lent to finite transformation (|3|). Functional equation (0) follows from the characteristic 
equation for the operator R, and for a given function f3, an explicit expression for G can 
be constructed by solving the corresponding Lie equations for operator (£|), 

- d9> dX (7) 



with the boundary conditions £'\\=o = £ and g'\x=o = 9- 

It is important in what follows that the invariance of the function / w.r.t. the RG 
transformations written as Eq. (|5|) is equivalent to the vanishing condition for the coor- 
dinate k of operator (^), which is written in the canonical form as 

R = aedj , se = I e - (3{g)I g = . (8) 

This condition must be considered on a particular solution I(£, g) of the boundary- value 
problem. 

The "exponentiated" variables 

x = e\ a = e\ g(a, g) = G(X, g), 

are natural in QFT. In these variables, the group transformation T(A) = T a becomes 

T a : {x' = xja , cf = g(a, g) } , g(l,g) =g, (9) 

and 

9(x,g) = g(x/a,g(a,g)) , T a ■ T b = T ab . (10) 

Functional equation ([HI and transformation (^) appear, for example, in a QFT with one 
coupling constant in the massless limit, where the dimensionless constant x = Q 2 / fi 2 is the 
ratio of the squared transferred four-momentum Q to the squared "normalized" momen- 
tum fi, g is the coupling constant, and the invariant g is the invariant (or effective) coupling 
function. More complicated QFT RG transformations are generalizations of transforma- 
tion (0) obtained by "multiplying" solution characteristics, g — > {g} = (gi, g 2 , . . . , g k ), 
and by introducing additional parameters in the functions (ji(x, {g}) (see Sec. 49 in 
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In the QFT model with two coupling constants gi = g and g 2 = h and one particle 
mass m, for example, the transformation T a becomes 

T a :{x' = x/a, y' = y/a, cf = g(a,y; g, h) , ti = h(a,y; g, h) } , (11) 

g(l,y;g,h) = g , h(l,y;g,h) = h, y = Q 2 / 'u 2 . 

In the case of several coupling constants, the characteristic equations that are connected 
with the infinitesimal operator 



R = xd x + yd y - Px(y; g, h)d g - (3 2 {y\ g, h)d h 



(12) 



Pi{y,g,h) 



9g(£,y,g,h) 



p2(y;g,h) 



dh(£,y;g,h) 



become the system of first-order PDEs 
dg(x,y;g,h) y _ 



x ■ 



dx 



,, _ dh(x,y;g,h) y _ - 

m-;g,h), x — = /3 2 {-;g,h) , (13) 



however, all these one-parameter transformations are based on the transformation of a 
single independent argument x. 



2.2 Self-similarity and the FS 

In a particular case where the function g is linear in its second argument, G ~ g, solutions 
of Eq. fllOD have a powerlike dependence on the argument x, i.e., g(x,g) = gx k , where k 
is a number, and transformations @ become transformations of powerlike self-similarity 
(the so-called automodel transformations) 

P a : {x' = x/a, g' = ga k ) , 

which are commonly used in problems in gas and liquid dynamics. From this standpoint, 
transformations ([5]) and ( |T0D for arbitrary G, g are functional generalizations gx k — > g(x, g) 
of customary self-similarity transformations and can therefore be called the functional 
self-similarity transformations ; this term is a synonym for the RG transformations. 

Therefore, in the RG-transformation framework, FS reflects the group nature of func- 
tional relations. The universality of the FS formulation in the QFT and classical physics 
models is due to the common scaling transformation and common functional transforma- 
tion of the solution characteristic = g(fi,g). Various realizations of the RG differ only 
in the form of the function (3(g), which is customarily calculated using an approximate 
solution obtained, as a rule, within the perturbation theory (PT). Therefore, the main 
role of the FS until recently was to establish a posteriori that a system under consider- 
ation admits functional transformations with a group structure. Group transformations 
themselves followed from additional considerations concerning system solutions [§]. 
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2.3 Constructing the RGS 

The situation changed substantially after passing to investigating boundary-value prob- 
lems in mathematical physics. The FS property of a desired solution here results in 
functional relations (which follow from the group nature of a solution) or in symmetries 
of the RG (=FS) type in the infinitesimal form. In contemporary mathematical physics 
models, RGSs can be regularly found by using the scheme in ||, which naturally incor- 
porates the FS and its formulation as the invariance condition for a particular solution of 
a boundary-value problem. 

To demonstrate the role of FS in constructing RGSs, we recall the four main steps 
in this construction. The first step is to construct a special RG manifold (differential, 
integral-differential, etc.), which differs, generally speaking, from the manifold determined 
by the initial system of equations describing the physical system under investigation. The 
second step is to find the maximum extended transformation group Q admitted by the 
RG manifold. The third step is to restrict the obtained group Q on a solution (exact 
or approximate) of a boundary-value problem. The transformation group that appears, 
which is just called the RG, is a set of infinitesimal operators Ri, each of which contains 
a solution of the boundary-value problem in its invariant manifold. The fourth and final 
step is to use the infinitesimal RG operators Ri to construct finite transformations of 
the group and to obtain an analytic expression for the solution of the boundary-value 
problem. 

The first and principal step in this scheme for RGS construction can be realized in dif- 
ferent ways depending on the mathematical model structure and on the type of boundary 
conditions || . The desired RG manifold can be obtained by adding parameters that enter 
the initial equations and boundary data to extend the list of group variables or appending 
to this list derivatives of the given dependent variables or nonlocal variables. We can also 
increase the number of initial equations by writing boundary data either in the embedding 
equation form or as additional differential constraints. Sometimes, the admitted group 
can be extended by dropping small parameters to simplify the initial equations. 

The maximum group Q, which is calculated for the RG manifold using the modern 



group analysis algorithms |T0[ in the second step, is not yet the RG, because it is not 
connected with a particular solution. 

The third step constructs the RG itself; namely, we restrict the obtained group Q 
on a particular exact or approximate solution of the selected boundary-value problem. 
Mathematically, the restricting procedure implies "combining" the canonical coordinates 
of operators of the group admitted by the RG manifold. The vanishing condition for the 
sum of these coordinates on the solution of the boundary-value problem (the condition 
of invariance w.r.t. the RGS operator, which is analogous to condition (|8|))) results in a 
system of algebraic identities, which relate coordinates of different operators and therefore 
generate the desired RGSs. From the physical standpoint, the verification of this condition 
in each actual case is the realization of the FS principle, which therefore establishes 
the invariance of the particular solution w.r.t. the RG transformations. On the other 
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hand, this step substantially relies on the PT solution of the problem under investigation 
simultaneously realizing the general principle of the RG method consisting in constructing 
an improved (in comparison with the PT) solution. In this scheme, the PT parameter 
can be any physically appropriate parameter or a set of such parameters. 

Verification of the FS condition in the RGS construction scheme above differs from the 
procedure for verifying the invariance of the analyzed solutions w.r.t. an iterated sequence 
of scaling transformations, which are used to analyze the behavior of a physical system 
in the theory of critical phenomena || and are also called RG transformations although 
they do not constitute a group in the general case (in contrast to the QFT models, where 
the group property is especially proved J?], §[]). This scheme (we call it the Wilson RG) 
is now widely used in mathematical physics to analyze the asymptotic behavior of DE 
solutions JTT], [T2j and to construct the envelope of the solution family ]T3| . When applying 



the Wilson RG to mathematical physics problems, the algorithm for improving the PT 
solutions that contain singularities consists of introducing additional parameters in the 
solutions, using these parameters to remove singularities, and demanding the solutions 



to be independent of the way these parameters are introduced [11]. Whether such a 
construction is consistent with the transformation group of the desired boundary-value 
problem solution remains an open question although this algorithm works successfully in 
some particular cases. 

The RGS in mathematical physics is defined by a set of RGS generators. Therefore, 
as in the classical group analysis of DEs, it suffices to consider only infinitesimal trans- 
formations that are characterized by an infinitesimal operator algebra. We note that this 
algebra often consists of more than one operator in contrast to the QFT-type models 
where we typically have only one RG operator. Both the dimension and the construc- 
tion of the RG operator algebra depend on the mathematical model and on the type of 
boundary conditions. 

The RGS group, which is restricted on a solution of a boundary-value problem, can be 
not only a point Lie group but also a Lie-Backhand group, an approximate transformation 
group, a nonlocal symmetry group, a non-Lie symmetry group, etc. ||. We present the FS 
conditions and various forms of the FS transformations pertaining to concrete problems 
in which the RGS is used. 



3 The FS analysis of systems that are close to quan- 
tum field systems 

We begin with the FS transformations for systems described by ODEs and, for simplicity, 
consider a boundary- value problem for the function u(t) satisfying a first-order ODE [|T^, 
14 , § with the parameters b k (k = 1, 2, . . .) explicitly entering the equation: 

u t = f(t,u,b k ), (14) 
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u{t) = x . (15) 

This example is not only methodologically important; it illustrates the FS application to 
QFT where such equations appear in the RG method when analyzing a system of DEs 
for invariant coupling functions. 

Constructing the RGS for boundary- value problem (0), ( [Top using the algorithm 



in m consists in adding to (|14D the embedding equation, which has the form of a linear 
first-order PDE, 

u T + f(T,x,b k )u x = 0. (16) 

The system of equations (|14|) and ( |16"D determines the desired RG manifold in the space 
of all group parameters {u, t, r, x, b h }, which are not only dependent and independent 
variables and the parameters of initial equation (PD but also the parameters r and x 
entering the boundary data. 

The symmetry group Q admitted by manifold (0), ( |T6"D depends on the form of the 
function /. In the practically important case (the ultraviolet limit of QFT models) where 
this function does not depend explicitly on time t and only three parameters b 1 = a, 
b 2 = b, and b 3 = c enter the initial equation, i.e., where / = f(u,a,b,c), the admitted 
group Q is determined by the seven-term operator 

7 

X = ^a,X i , (17) 

i=l 

The first two functions a,\ and «2 in (|17]) are arbitrary functions of all seven group variables 
{t, t, x, a, b, c, u}, and the remaining functions depend arbitrarily on the parameters a, b, 
and c and on the combinations t = t — {1/ f(u)) and f = r — (1/ f(x)). The angle brackets 
denote integrals w.r.t. the respective variable u or x. Explicit expressions for three of the 
seven operators entering ( |17D are 



X 1 = d t + f(u)d u , X 3 = f(u)d u , 

X, = f{x) < f a (x)/f 2 (x) >d x + f{u) < f a {u)/f\u) > d u + d a . 

The remaining operators are obtained by substitutions: X 2 and X 4 are obtained from X\ 
and X% with the respective substitutions t —> t and u — > x, and X§ and Xj from X5 with 
the respective derivative substitutions d a — > db and d a — > d c . 

The operation of restricting the group Q is the verification of the imposed FS condition, 
which is analogous to the equality x = in (|8|), i.e., the solution must be invariant w.r.t. 
the RG transformations, or, in other words, the coordinate of the canonical operator X 
must vanish on the solution of the initial problem u = U(t,x,r,a,b,c). If a restriction 
of the group Q admitted by manifold (|14D , ([16]) on an (approximate) solution of problem 
(PI), flT5| ) is fulfilled, the "restricted symmetries," which we call RGSsJ^| appear. 

2 Some R operators of RGSs thus denned may result from symmetries of the equations. 
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Which perturbation series becomes the function U depends on the actual problem 
setting. For instance, for a polynomial function /, 



f = au 2 + bu 3 + cu 4 (18) 

we can choose the PT over the variables (t — r) or over the parameter a, b, or c. Examples 
of restricting the group on a solution with a PT in a with b = c = was considered 
14| . ||, and with a PT in b with a = 1 and c = in 0. We present two possible RGS 



in 



operators appearing after such restrictions JTJ 



R x = x 2 rd x + d a + u 2 td u , b = c = 0, 

(19) 

R 2 = (a; 2 (l + bx)r + x) d x + (m 2 (1 + 6tt)t + u) <9„ - bd b , a = 1 , c = . 

We can now obtain the solution of the Cauchy problem if we write the condition for its 
invariance w.r.t. the RGS operator, namely, the FS condition, which is a linear first-order 
PDE. For instance, for the operator Ri, we have 

tu 2 - x 2 tu x -u a = 0. (20) 

Solving the characteristic equations (the Lie equations) for this equation, we obtain the 
desired solution of the boundary- value problem for b = c = 0. The FS condition of 
type fl20D is used twice: first, when constructing the RGS operators, we substitute an 
approximate solution U for u in (|20|) ; second, when finding the solution of the boundary- 
value problem, we use these RGS operators. 

In both examples, constructing the PT is easy and results either in the powerlike 
dependence or in a combination of powerlike, reciprocal, and logarithmic dependencies 
on the initial data. In this algorithm, constructing the PT in powers of b, which gives 
the second operator R2 in flT9]), starts from the unperturbed state, which is chosen to 
be the solution of the boundary-value problem at a = 1 obtained by applying the first 
operator Rx, i.e., in the approximation of strong nonlinearity (w.r.t. the parameter a). 
Therefore, the improvement of the PT using the operator R 2 is reduced to the consequent 
improvement of the PT for the boundary- value problem with the function / = au 2 + bu 3 
first in the parameter a and then in the parameter b using the respective one-parameter 
groups R\ and R 2 . 

In many cases, such a procedure of improving the PT consecutively over several pa- 
rameters of a model fails. One possible reason is the absence of a "nonlinear" unperturbed 
solution. For the function / discussed here, which is chosen as a polynomial in u, con- 
structing the PT in c with fixed a and b is not so simple as constructing the PT at c = 0, 
because in the case with fixed a and b, the unperturbed state is determined in terms of 
the Lambert function |[6| , which admits no simple analytic representation (see also [17) ). 



Therefore, constructing an RGS that generates group transformations w.r.t. several 
parameters simultaneously and improves a PT that admits an easily found unperturbed 
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solution is important. In fact, we want to construct a two-parameter RG. This problem 
will be discussed in detail elsewhere; we consider only one example of such an RG here: 



R3 = -< h{u)/f{u) > d t - < f b (x)/f\x) >d T + d b , 

(21) 

i?4 = " < fc(u)/f(u) > d t - < f c (x)/f(x) >d T + d c . 

Verifying the FS conditions, i.e., restricting the group Q on the PT solution over the two 
parameters b —>■ and c — > simultaneously is not difficult technically, and using the finite 
transformations generated by the operators R3 and R4, which generate a two-dimensional 
algebra, we obtain the desired solution of problem fll4|) , Ql5[) . Therefore, using the two- 
parameter RG, we can avoid mathematical problems that arise when using one-parameter 
RGs. 

We conclude this section with a few remarks. We have demonstrated the use of the 
FS condition to construct RGSs that improve the PT over the parameters a, b, and c 
entering the equation, although the corresponding operators by no means exhaust an 
infinite set of RGSs, which is parameterized by a continuum set of the RG operator 
coordinates |IR fiM M. Quite analogously, RGS operators that improve a PT w.r.t. an 



appropriate combination of dynamic variables and initial data (e.g., w.r.t. the difference 
t — t) can be found. Examples of such RGSs are presented in the following sections. We 
only mention here that the FS condition in the form of equality (|8|) is obviously analogous 
to FS condition ( p0|) pertaining to boundary- value problem (|14D , (|I5D . The RGS examples 
and the corresponding FS condition considered in this section follow from a single ODE. 
Similar constructions are valid for systems of ODEs depending on several parameters. 
Examples of embedding equations for such systems, which were used to find RGSs, can 
be found in [14]. 



4 The FS conditions for PDEs 

In this section, we consider systems that are described by mathematical models based on 
PDEs or systems of PDEs. In contrast to the previous section, where the appearance of 
first-order partial derivatives in FS conditions is due to a transformation of parameters 
in initial equations and due to embedding equations, which imply taking the boundary 
conditions into account, we here demonstrate boundary- value problems for which the FS 
conditions are systems of DEs (constraints) that contain higher-order partial derivatives 
w.r.t. the independent variables. The differential formulation of the FS conditions then 
permits constructing the desired solution of the problem; finding finite FS transformations 
results in formal power series (see fT8| , Vol. 3, Chap. 1). In addition, we give an example of 
an FS condition (simpler than in Sec. 3) in the ODE form whose solution also reconstructs 
the solution of the boundary- value problem from the given PT. 
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4.1 Boundary- value problem for the Burgers equation 

An explicit example, which illustrates the variety of the FS condition formulations for 
systems based on PDEs, is provided by the boundary-value problem for the modified 
Burgers equation, 

u t - au 2 x - vu xx = , (22) 

u(0,x) = f(x), (23) 

with the nonlinearity parameter a and dissipation parameter v. The continuous point 
symmetry group admitted by manifold is defined by nine operators. Six of them 



are symmetries of the equation and have been discussed in the literature (see, e.g., |18 
Vol. 1, Chap. 1, p. 183). They correspond to the projective transformation, the dilation 
transformation in the plane (t, x), translations along the axes t, x, and u, and the Galilean 
transformations. The seventh operator is the operator of an infinite Abelian ideal = 
a exp(— au/v)d u of the group; the coordinate of this ideal is parameterized by the function 
a(t, x, a, v) of four group variables restricted by the linear parabolic equation 

at - va xx = 0, 

which coincides with the linear part (a = 0) of initial equation (0). Eventually, when we 
interpret the equation parameters as independent variables, two more operators appear 
and involve these parameters in the group transformations; these operators correspond to 
scaling transformations of the respective variables a and v [|l9|j . 

Restricting the group admitted by manifold ( p2|) on the solution u = U(t, x, a, v) of the 
Cauchy problem, i.e., verifying the FS condition, results in an algebraic relation, which 
expresses the coordinate of the infinite-dimensional subgroup generator (the function a) 
through the coordinates of the remaining eight operators at any time t, including t = 
when this solution U (0, x, a, v) = f(x) is known from boundary condition Q2"3|). Using the 
standard representation for the solution of a linear equation on the function a with the 
initial value a(0, x, a, v) obtained from the FS condition and substituting this representa- 
tion in the formula that determines the general element of the Lie algebra, we obtain the 



desired RGS operators. Thus, we obtain the RGSs for boundary- value problem (f22|), (|23|) 
by combining the symmetries of the eight-dimensional algebra generated by the above 
operators and symmetries of the infinite-dimensional subalgebra generated by the opera- 
tor Xqo. Each of the eight RG operators obtained (and their linear combinations whose 
coefficients are arbitrary functions of a and v) contains a solution of the Cauchy problem 
u = U(t, x, a, v) in the invariant manifold and permits finding group transformations of 
both the variables {t, x, a, u} and various functionals (local and nonlocal) of the solution. 

We present two such RGS operators, which improve the corresponding PTs over the 
parameter a (operator R 5 ) and the independent variable t (operator R G ), 

Rs = d a + - a (-u + exp (-^) < f(x) >) d u , (24) 
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i?6 = d t + exp(-^) <afl + vf xx >d u . (25) 

The double angle brackets here are integral convolutions of the corresponding functions 
with the fundamental solution G(t, x, v) of the linear equation for the function a multiplied 
by the exponent of the function / from boundary condition (|23|): 



— oo 

The invariance conditions for a solution of a boundary-value problem (the FS conditions) 
corresponding to RG operators (EMI) and are the two first-order ODEs 



-u a -- + -exp( )</(a;)>=0, 26 

a a \ v J 



-u t + exp(-—^<af* + vf xx >=0. (27) 

Approximate solutions of Eq. (p2p can be extended in parameters a or t by solving either 
Eq. (p6|) or (1271) with the proper initial conditions. This eventually leads to the common 



exact solution u = {y/a) ln((l)), which is valid for all values of the parameters a and t flT5 |. 

In the example above, FS condition (|26| ) for the boundary-value problem for a PDE 
becomes a first-order ODE. Formally, this equation is simpler than the FS conditions for- 
mulated for the ODE boundary-value problems in Sec. 3 and differs from the corresponding 
FS conditions for QFT models. On the other hand, formulating boundary conditions for 
a first-order ODE in the embedding equations language, we notice that condition (EJ 



can be treated as an embedding equation for boundary- value problem (p22|) , (|23|) with 
the embedding parameter a. In turn, the invariant embedding method can also be ap- 
plied to boundary-value problem (|26|) with the variables {a, u} and the initial condition 
u = uo(t,x, v) at a — 0. Then, because the initial value uq depends on the "parameters" 
{t,x,v} entering ODE (|26|) , the corresponding embedding equation becomes integral- 



differential JT51. The set of the RGS operators for boundary- value problem (EZ2]), ( |2"B| ) 
includes not only the above FS conditions in the form of first-order ODEs but also opera- 
tors that result in FS conditions both in the form of first-order PDEs, which are analogous 



to the operators previously considered, and in the form of mere algebraic relations [{19 
Therefore, various FS can relate to the solution of the same boundary-value problem. 



4.2 The boundary-value problem for nonlinear optic equations 

It is sometimes impossible to make a boundary-value problem PT and the FS conditions 
consistent if we confine ourselves to only point symmetries. Below, we consider an example 
where the FS condition is a second-order PDE. We now consider the following boundary- 
value problem for the system of equations for the nonlinear optics of collimated wave 
beams: 

v t + vv x — an x = , n t + vn x + nv x = , (28) 
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u(0,a;)=0, n(0,x)=N(x). (29) 

Here, the dimensionless coordinates t and x describe the spatial evolutions of the deriva- 
tive v of the beam eikonal and of the dimensionless intensity n in the direction into the 
bulk and in the transverse direction, and a is the parameter of the nonlinear refrac- 
tion. The hodograph transformation reduces (|28|) to the system of linear equations in the 
variables t = nt and x = x — vt 

r w ~ n\ n = , Xw + ocT n = , (30) 

r(0,n) = 0, X (0,n)=H(n). (31) 

where we use the notation w = v/a. Formal construction of the RGS and analysis of the 
FS conditions for system ( |3Tj ) could be as for the Burgers equation. As previously, a finite- 
dimensional subalgebra of point symmetry operators and its infinite-dimensional ideal 



arise. However, the procedure of restricting this ideal in order to obtain the RGS |20|, £1 
implies solving the system of linear PDEs, which coincides with initial system (|30D; this 
is an intrinsic feature of linear equations. 

Therefore, we construct the RGS using not the point symmetries but the Lie-Backlund 



symmetries (see |[L8| , Vol. 3, Chap. 1), for which the terms "higher" and "generalized" sym- 
metries are also used. The desired RGS operators for boundary- value problem (|30|), (|3l| ) 
can then be conveniently written in the canonical form 

R = fd T + gd x , (32) 

where the coordinates / and g are functions in the extended space of group variables, 
which, in addition to the set (r, %, w, n, a), includes derivatives of r and x °f arbitrary 
finite order in n (and, perhaps, in a). The problem of constructing the Lie-Backlund 
symmetries for boundary- value problem (|30"1), ( 3l|) was discussed in detail in | 20|| . It was 



shown that these symmetries are generated by operators of form (|32| ) whose coordinates 
are linear combinations of r, x-, an d their first- and second-order derivatives depending 
on w and n. We write the expressions for the coordinates of only three operators: 

fi = -r/2 + nr n + (l/2)nwxn, 9i = -(aw/2)r n + n%„; 



f2 = nr n , g2 = Xn + nXnn, 

(33) 

/ 3 = (l/4)r - nr n - (5/4)w«x„ + (-n + (a/A)w 2 ) nr nn - wn 2 Xnn, 

g 3 = (3/4)fT„ - (2n - (a/A)w 2 ) Xn + awnr nn + (-n + (a/A)w 2 ) nx nn ■ 

The first operator with the coordinates fi and g\, which are linear in the first derivatives, is 
equivalent to the point symmetry operator; the other two operators are the Lie-Backlund 
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symmetry operators of the second order. Restricting the Lie-Backlund group found for 
the RG manifold consists in verifying the FS condition 

/ = 0, = 0, (34) 

which must be satisfied on the solution of the boundary-value problem and be consistent 
with boundary conditions fl30l) . The form of the coordinates / and g implies that these 
conditions, written in terms of the boundary function H(n), must satisfy linear PDEs 
with variable coefficients. In particular, for the case N(x) = cosh" (x), the FS condition 
is fulfilled for the combinations f = fi + 2(/ 2 + ^3) and g = gi + 2(g 2 + #3), and the desired 
Lie-Backlund RGS operator is 

R 7 = (2n(l - n)r nn - nr n - 2nw(x n + nx nn ) + t| nw 2 T n j d T 

(35) 

+ (2n(l - n)xnn + (2 - Sn)x n + aw (2nr nn + r n + ^ (nXnn + Xn))) d x 



Condition fl3~4] ) for the operator R? results in a second-order PDE, which distinguishes 
it qualitatively from the first-order equations arising in the FS conditions considered 
previously. Namely, the finite FS transformations for the latter can be found in a closed 
form as the solutions of the corresponding Lie equations expressing the transformations of 
dependent and independent variables and parameters. The solutions of the Lie equations 
for the RGS with Lie-Backlund operator (^) result in finite transformations written as 
formal power series. However, this does not mean that FS conditions have no practical 
importance in this case; on the contrary, the differential formulation of the FS conditions 



permits considering relations (|34T ) as additional differential constraints, which must be 
taken into account when finding solutions of the boundary-value problem. The solution 
of the initial equations with the differential constraints taken into account leads to the 



solution of the boundary-value problem j2(J, ^lj, ^] that is invariant w.r.t. RGS group (B5| 



To conclude this section, we note that the FS conditions written as vanishing conditions 
for coordinates of a Lie-Backlund canonical RG operator, being considered as a system of 
differential constraints, play the role of embedding equations. Imposing the FS conditions, 
we come to the problem of finding the solution that is invariant w.r.t. the found RGS 
operator. A more informative example is provided by using the FS conditions to find 



point RGS operators as the group admitted by the initial equations and FS relations [JT5 
There, RGS operators containing higher-order derivatives w.r.t. parameters entering the 
equations can be discussed in principle. 



5 The FS and RGS for systems with small parame- 
ters 

We now consider the FS property for systems described by DEs with small parameters. If 
a system contains a small parameter a, we can start the RGS construction by considering 
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the simplified (a = 0) model, which admits a wider symmetry group in comparison 
with the case a / 0. When we take the contributions from small a into account, this 
symmetry is inherited by the initial system of equations, which results in the appearance 
of additional terms, corrections in powers of a, in the operator coordinates. Restricting 
the obtained symmetry group to an exact or approximate solution of the boundary-value 
problem, we obtain the desired RGS, which can be also represented as operators whose 
coordinates are infinite series in powers of small parameters. For boundary data of a 
special form, these series terminate, which results in an exact RGS without restrictions 
on values of the relevant parameters. For boundary data of a general type, the procedure 
of truncating the infinite series to a finite number of terms leads to an approximate RGS 
at small parameter values. 



5.1 FS conditions in nonlinear plasma theory 

We now consider the boundary- value problem for the equations of the nonlinear interaction 
of laser radiation with a plasma |^v|] . Such an interaction for a p-polarized electromagnetic 
wave of frequency u, for which only the z component of the magnetic field is nonzero and 
which propagates from the vacuum toward an inhomogeneous plasma, is described by 
a system (we do not present it here) of nonlinear nonstationary equations for six scalar 
functions: two components of the electron velocity v x and v y , the electron density n, two 
electric field components E x and E y , and the z component B z of the magnetic induction; 
these functions depend on the time t and the two coordinates x and y. 

The nonlinearity of these equations is essential in a small space domain near the 
plasma resonance (at uj\ ~ u 2 ) where the presence of natural small parameters (such 
as the smooth inhomogeneity of the ion density N(x) along the x axis and the small 
incidence angles $ of laser beams at the plasma) implies the appearance of a hierarchy of 
components of the p-polarized light wave at the critical plasma point. When constructing 
the inherited point RGS, this allows reducing the total system of six initial equations to a 
simpler system of two one- dimensional nonlinear PDEs for the components of the electron 



velocity v and the electric field E along the density gradient vector p4 , [23 



2 2 47T \eei \N 

v t + avv x — E = 0; E t + avE x + uj l v = , to L = . (36) 

m 

Here, the functions v and E are expressed in units of the dimensionless nonlinearity 
parameter a, which is proportional to the value of the magnetic induction B at the 
critical point at the laser frequency; the coordinate y enters only in combination with 
time, t — > t — ($y/c). 

The system of six initial equations admits only a finite group of point transformations, 
namely, the group of translations along the t and y axes for arbitrary N(x). At a constant 
ion density, N = const, the additional group of a:-axis translations and the group of 
simultaneous rotations in three planes, which are determined by the coordinates {x, y} 
and the corresponding x and y components of the electron velocity and of the electric 
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field, arise. In contrast to the initial equations, system ( PB"| ) admits an infinite group of 
point transformations with the operator containing the three terms 

X x = mY ; X 2 = ^ 2 d x + - Y(fi 2 )d v + - Y\fi 2 )d E ; 

a a 

(37) 

^3 = - A*3 ( a da - vd v - Ed E ) ; Y = d t + avd x + Ed v - uj 2 L vd E , 
a 

each of which contains an arbitrary function /ij of independent and dependent variables 
and of the parameter a while the differential constraints 

y 3 W + y(4) = o; y(^) = o. 

are imposed on \i 2 and /i 3 . 

To obtain the RGS, point transformation group ([37]) must be restricted on a solution of 
the boundary-value problem that is approximate over the powers of the parameter a; this 
solution is such that the leading approximation for the functions v and E is determined by 
a solution of the linearized system of initial equations with the corresponding boundary 
conditions (the propagation of the electromagnetic wave from the vacuum toward the 
plasma) and with the given shape of the density N(x) in the plasma resonance domain 
taken into account; the corrections proportional to a appear when linearizing system (|36|). 
Verifying the FS conditions for group fl37D , we find that fii = 0, fi 2 = —E/u 2 , and fi 3 = 1 
for this particular solution, which gives the desired RGS operator 

R B = X 2 + X 3 = ~d x + d a . (38) 
us 1 

The condition of invariance of the solution of the boundary-value problem w.r.t. the 
RG operator R s (the FS condition) is described by the system of first-order PDEs 

v a -^v x = 0, E a -^E x = 0. (39) 

The solution of the characteristic equations for this system (the system of Lie equations 
for the operator R$) reconstructs the PT in the parameter a up to the exact solution of 
the boundary- value problem p3[, 



E = — [qi sin ojt + q 2 cos cut) , v — — \q\ cos uit — q 2 sin ut) , 

A A (40) 

x = fj, + e(qi sin cut + q 2 coscut) , e = aL 2 / A 2 , 

which describes the nonlinear structure of the electric field in the plasma resonance do- 
main. The dimensionless quantities x and /x are normalized to the resonance width A, and 
L is the characteristic inhomogeneity scale of the plasma ion density. The choice of the 
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form of the functions q\ and q 2 in relations (|4T)| ) is determined by the actual dependence 
on the x coordinate of the electric field, which is a solution of the initial total system of 
equations under the corresponding boundary conditions and with the given shape of the 
density in the plasma resonance domain. 

For example, for a cold electron plasma with a linear density function, the functions q\ 
and q 2 are 

g 1 = (l+ /U 2 )- 1 , g 2 = /i(l + /i 2 )- 1 , (41) 

where the width A = [yjuS)L is determined by the frequency of plasma particle collisions. 

For a hot plasma in which the heat motion of electrons is essential, linearizing the 
initial system of equations, we find that in the case where the density function is linear 
in the plasma resonance domain, the electric field distribution is expressed through the 
Airy-Fock functions. The corresponding nonlinear structures of the electric field and 
of the density and velocity of electrons, which appear from using RG operator (p8|) to 
continue the linear relations with the first nontrivial heat correction (for a — > 0) to finite 
values of a, are determined, as before, by expressions (f40|) in which the resonance width A 
now depends on the electron temperature, A = (^V^L/co 2 ) 1 ^ 3 , and the functions qi are 

oo oo 

q x = J dti cos(/i£ + £ 3 /3) , Q2 = J sin(^ + £ 3 /3) • (42) 



A feature of the formulas for v and E in ([f(]) with the functions q\ and q 2 from fl42|) is 
that they give the exact (at uj\ = u 2 ) solution of Eqs. (|36D , in which the electron pressure 
is neglected but the nonzero electron temperature is nevertheless taken into account. 
These formulas as well as their physical consequences were analyzed in detail in [25 



Here, we only note the value of the results obtained: the FS condition for a simple mathe- 
matical model leads to proper results even in the leading order in which small corrections 
to RGSs are neglected (although a modification of RG operator ( |5BD taking these correc- 
tions into account is not complicated; the corresponding expression for corrections in the 
density gradient was given in f23|). The big freedom in choosing the functions q\ and q 2 



in Eqs. ([!(]) permits analyzing the nonlinear structure of the electric field in the plasma 
resonance domain in both cold and hot plasmas uniformly. 

We note that RG operator (|38| ) is qualitatively analogous to Bogoliubov RG oper- 
ator (|j). As in the standard QFT RG case, operator (^) determines the translation 
transformation along a solution characteristic (the parameter a) as well as more compli- 
cated functional transformations of the coordinate x. The FS property of the physical 
system under investigation is the invariance of the desired functions, the field E and elec- 
tron velocity v, w.r.t. these RG transformations. Then, the nonlinearity parameter a, 
which enters the RG transformation, is not assumed to be small in contrast to the terms 
dropped when obtaining the operator R$ of the "approximate" RGS; the smallness of 
these terms was crucial for passing from the initial system of equations to Eqs. (p6|) . 
Therefore, the RGS can be used to extend the solution obtained in the form of a power 
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series in the parameter a to the exact, essentially nonlinear solution, Eqs. ([5])- Then, as 
in the QFT case, the obtained RG operator corresponds to the exact group transformation 
w.r.t. the parameter a for a solution of Eqs. (|36|). At the same time, this RG operator is 
approximate because it was obtained by neglecting parameters other than a and admits 
a subsequent improvement w.r.t. those parameters. 



5.2 The FS in nonlinear wave optics 

In many cases, we cannot neglect all but leading (zero-order in small parameters) contri- 
butions to an RGS. Therefore, we consider a variant of the RGS and FS conditions that 
explicitly takes contributions of small parameters into account in the RG operator; such a 
modification appears when analyzing the boundary- value problem, which generalizes Q2"5|) 
and (||), __ 

v t + vv x - an x - (3d x {{x l ~ u / ' ^fn) d x (x u ~ x d x [\fn))) = , 

(43) 

n t + nv x + vn x + iv — 1) — = . 

x 

Here, a is the nonlinear refraction parameter as previously, (3 is the parameter that defines 
the diffraction effects, and the respective values v = 1 and v = 2 correspond to flat and 
cylindrical geometries of the beam. Boundary conditions for Eqs. (|43| ) determine the 
curvature of the wave front of the beam and its intensity distribution over the transverse 
coordinate x, 

v(0,x) = V(x) = -x/T, n(0,x)=N(x). (44) 



5.2.1 Flat geometry 

We now set the geometry to be flat (u = 1), neglect the diffraction ((3 = 0), and consider 
the simplest case where the wave beam on the medium boundary (i.e., at t = 0) has a 
flat wave front. Then, Eqs. (f£5Q become system We write the RGS operators for 

boundary- value problem (|30|) , (|3i|) in canonical form (|32]) . We represent the coordinates / 
and g as power series, 



/ = 5>T; = X>V- (45) 



i=0 



The method for calculating the coefficients f l and g l and the resulting system of recursive 
relations for these coefficients was derived in [^, ^7|. There, each coefficient was deter- 



mined up to an arbitrary function of n, X(s)> an d of the combination 7"( a )— w(sx( s )+nx{s+i)) 
where (s) denotes the sth-order derivative in n. This arbitrariness can be removed by 
restricting the group on the desired solution of the boundary- value problem, i.e., by im- 
posing FS condition (|34|), which becomes a differential or algebraic constraint that is 
compatible with boundary condition ([U]) at r = 0. Truncating series (|45|) , which is 
possible at small values of a, we obtain an approximate symmetry; the contributions f° 
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and g°, which do not depend on the parameter a, are then determined by the system of 
equations 

r w - nxn = , Xw = . (46) 

which is simpler than initial system fl5U|). In contrast to system ([H]), which admits only 
a finite group of Lie-Backlund transformations of a given order, system fl46|) admits an 
infinite symmetry group because the coordinates f° and g° can be arbitrary functions of 
their arguments. At small a, the symmetry of Eqs. ( [46] ) is inherited by initial system (|30| ) 
up to an arbitrary given order in a. Restricting the obtained approximate group on the 
solution of the boundary- value problem, we obtain the desired RGS. 

We now turn to RGS operators with coordinates in the form of the binomials 

f = f° + af\ g = g + ag K (47) 

These operators depend on the functions f° and g°, which, in turn, can be (nonuniquely) 
determined by the function N(x) of the transverse intensity distribution at the boundary 
of the nonlinear medium. 

For the "soliton" beam intensity distribution function N(x) = cosh - x, we can ob- 
tain two sets of formulas for the coefficients resulting from different expressions for f° 
and g° [|(| : 



a) f° = 2n(l - n)r nn - nr n - 2nw(x n + nXnn) , f 1 = - nw 2 r nn , 

g = 2n(l - n)xnn + (2 - 3n)x„ , g 1 = w {2nr nn + r n ) + — (nx nn + Xr. 



(Ah 



b) f° = 1 + 2nxn tanh x , f 1 = (— ~ 2rr„ + 2r 2 tanh x] cosh -2 x , 

\n J (49) 



g° = , g 1 = —2rxn cosh 2 x ~ 2r n tanh x ■ 

Calculating contributions of higher orders |26|, |2^, ^| , we find that the functions f l and 
g % vanish for i > 2 in case a and formulas fllffi) and ( ^5|) describe an exact RGS; the 
comparison with coordinates (|35| ) demonstrates that the functions f 1 ' 2 and g 1,2 are the 
coefficients of the coordinate expansions in powers of a. In case b, series ( ^5|) do not 
terminate, and formulas (|49| ) and ( 45[) pertain to an approximate RGS. 

For the wave beam with a Gaussian initial intensity shape, N(x) = exp(— x 2 ), we have 
two sets of the / and g operators of an approximate RGS: 

a) f° = 1 + 2n X Xn , g° = , f 1 = -2rr n + — , / = -2 (r Xn + X r n ) , (50) 

n 

b) f° = 2n(rxn + r n x), g° = 1 + 2n X Xn , f l ^2 X r a , g 1 = 2 ( X Xa ~ rr n ) . (51) 
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Formulas (^)-(^TJ) demonstrate the main advantage of the approximate RGS method, 
which permits analyzing boundary-value problems with arbitrary boundary data, which 
are expressed through differential or algebraic expressions for the functions f° and g°. 
The example of operator (f48|) shows that in some cases the approximate RGS method 
can result in an exact RGS whose presence must be established using the tools in the 
previous section. Because expressions ( pE5|) for the coordinates f l and g l contain second- 
order derivatives, the corresponding operator R is the RGS Lie-Backlund operator of 
the second order, while operators (^9|)-(|51|) are equivalent to point symmetry operators. 
Meanwhile, because operator (|5l| ) contains the first derivative in the parameter a, the FS 
transformations also involve the nonlinear refraction parameter ||26|| . (See |26], ^] for a 



detailed analysis of formulas (^8|)-(pT|) and the physical consequences of their substitution 
in FS conditions 



5.2.2 Cylindrical case 



Analogously to the previous case, approximate RGSs for the case (3 ^ can be constructed 
using FS conditions [^7j. An example is the RGS operator for the cylindrical [y = 2) wave 
beam, 
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where the function S depends on the variable x = x ~ 

(3 



S( X ) = aN( X ) + 



Xx/Nix) 



contains two small parameters a and /3, and is determined by the initial beam intensity 
distribution function iV at the boundary of the nonlinear medium. Then, as in the case 
(3 = 0, there exist such functions iV for which operator ( |52] ) corresponds to an exact, rather 
than approximate, RGS, i.e., the RGS that is applicable at arbitrary, not necessarily small, 
values of the parameters a and f3 [27|. The FS transformations appear from the solution 
of the characteristic equations for the first-order PDE that is conjugate to (|52|) and permit 
continuing the PT solutions, which are determined only in a small vicinity of the nonlinear 



medium boundary, to a domain where essentially nonlinear effects prevail [[Z7 

To conclude this section, we again note that we can write the FS conditions in many 
different ways, which implies different approaches to their subsequent use in problems 
with a small parameter. The use of approximate symmetries results in some peculiarities 
in constructing the RGSs and in using the FS conditions. 
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First, as discussed in the previous section, we apply the FS conditions in the form of 
a system of higher-order DEs and analyze them together with the initial equations. In 



contrast to finite FS transformations, which are infinite formal series (see [p!q| , Vol. 3, 
Chap. 1), the differential formulation of the FS conditions is suitable for constructing 
solutions of a boundary-value problem. 

Second, the FS conditions expressed as approximate symmetries are additional differ- 
ential constraints, which together with the initial equations determine the RG manifold. 
In turn, such a manifold can be used to construct approximate point RGSs. 

Third, the FS conditions are important at all stages of constructing approximate RGSs 
already starting with setting the RGS operator coordinates for an unperturbed solution. 
New prospects are provided by the possibility of constructing RGSs based on approximate 
symmetries for problems with arbitrary boundary conditions. 

Fourth, we introduce parameters that are used to construct approximate RGSs in the 
set of variables of the FS transformations. 



6 Conclusion 

We now formulate how the FS notion has been transformed in the last decade. This 



notion first appeared as pertaining to the group transformations developed by Lie in ||29|| , 
which was used to construct the QFT (Bogoliubov) RG j|, The latter concept was 
based on the one-parameter Lie group of local transformations, the class the Bogoliubov 
RG belongs to. The FS notion unified various RGs and exhibited their intrinsic group 
structure. 

The progress in applying mathematical methods to the investigation of symmetries of 
differential and integral-differential equations, which resulted in modern group analysis 
(see, e.g., [lCj, [18], and references therein), has led to the transformation of the FS 

notion. 

First, this resulted in constructing the special class of symmetries related to the RG 
transformations, namely, the RGSs. Such a method (see Sec. 2) regularizes the procedure 
for finding RG transformations, at least for systems that are described by differential and 
integral-differential equations. Using the infinitesimal approach, we can formulate the RG 
invariance conditions through the RGS operators and can therefore use all the powerful 
tools of modern group analysis to construct these conditions. As a result, the RGS 
notion now includes not only point symmetries (Sec. 3) but also Lie-Backlund symmetries 
(higher, or generalized, symmetries) (Sees. 4 and 5), approximate symmetries (Sec. 5), 
etc. The list of examples, which is far from complete, presented in this paper includes 
the most advanced results in constructing the RGSs; furthermore, our approach is equally 
applicable to the cases of nonlocal RGSs and to RGSs of integral-differential equations if, 
for example, we apply the method described in |HJ. Therefore, modern group analysis 



is as important for developing the FS notion as was the classical Lie group analysis for 
establishing the Bogoliubov RG. 
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The FS conditions are formulated as conditions of invariance w.r.t. transformations 
determined by RG operators. This becomes important when the formulation of finite FS 
transformations is doubtful, e.g., for the Lie-Backlund RGSs (see Sees. 4 and 5). Being 
universal, this formulation is useful not only at the stage of finding the solutions but 
also when constructing RGS operators. Thus, the FS conditions and the procedure for 
their verification became important components of the RGS construction. The conditions 
themselves can be algebraic as well as differential relations containing derivatives of higher 
(not necessarily first, see Sec. 4.2) order. 

A feature of the new way of constructing RGSs is the regular procedure for constructing 
Lie algebras of finite dimension. Numerous examples of such algebras are given in H, 



19| , although they have not yet found wide practical applications (see Sec. 3). Such 
Lie algebras correspond to multiparameter groups. The RGS apparatus results in a set 
of group operators constituting a Lie algebra. Such groups are customarily used by 
researchers into symmetries of equations arising in contemporary theoretical physics. f\ 
However, the "renormalization group" in both the QFT (Bogoliubov) and the Wilson 
approaches was always assumed to be a one-parameter group. Rare attempts to consider 
two-parameter constructions have always (to the best of our knowledge JHJ) led to direct 
products of two one-parameter RGs. 

The RG technique was customarily used in problems with singularities to improve 
the PT and to give a correct description of the solution behavior in the vicinity of a 
singularity; these properties are intrinsic in our approach as well. Namely, using the FS 
conditions formulated in Sees. 4 and 5 on the base of the Lie-Backlund RGS, we were 
able not only to describe the structure of the known singular solutions but also to find 
p2| , |26| , p7j). We stress that the RGS approach to these problems results in 



new ones see 



a structure of two-dimensional singularities that differs from the structure of singularities 
appearing within the Bogoliubov RG setting. 
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